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OWTRAN is a FORTRAN computer program, developed at the U.S. Air 
Force Geophysics Laboratory (AFGL), which calculates the transmittance 
of the earth's atmosphere in the spectral region from 0.25 to 28.57um 
with 20 cm spectral resolution on a linear wavenumber scale. A choice 
of six atmospheric models covering seasonal and latitudinal variations 
from sea level to 100 km are available, in addition to a capability of 
allowing the user to input atmospheric profile or meteorological para- 
meters of his own choosing. The program includes four aerosol models 
and permits either hazy or clear haze condition to be selected, in addi- 
tion to the capability of selecting a particular visibility of the user's 
choosing. The model accounts for molecular absorption, molecular scat- 
tering, and aerosol extinction, plus atmospheric refraction and earth's 
curvature eer ee 

In examining AFGL's version of LOWTRAN (hereafter referred to 
AFGL/LOWIRAN), the author became aware that it would be possible to re- 
design the program without changing the basic methodology. This re- 
design achieves two goals: (a) a significant decrease in program size 
and hence computer time; and (b) a more structured program flow. The 
first goal is critical for repeated uses of the program or if the pro- 
gram is to be incorporated into a weapon/sensor system$ simulation. 
The second goal is tied to the first one in that a smooth, structured 
logic flow usually improves computer time. Moreover, a structured fro- 
gram is easier to understand, to use, and, when necessary, to modify 
and update. 


This report documents the updated code, in addition to including 
some corrections which have been brought to the author's attention. 
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FOREWORD 


This technical memorandum contains a brief description and program 
listing of the atmospheric transmittance computer code NWC/LOWTRAN. The 
computer code development has been supported by a contract with Naval 
Air Systems Command (AIR-503E). 


This document has been prepared primarily as an interim presentation 
of timely information. Although care has been taken in the preparation 
of the technical material presented, the results herein are to be con- 
sidered as preliminary in nature. 


P. B. HOMER 

Head, Weapons Systems Analysis Division 
Systems Development Department 

28 March 1977 
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INTRODUCTION 


LOWTRAN is a FORTRAN computer program, developed at the U.S. Air 
Force Geophysics Laboratory (AFGL), which calculates the transmittance 
of the earth's atmosphere in the spectral region from 0.25 to 28.57um 


(350 to 40,000 em”) with 20 can? spectral resolution on a linear wave- 
number scale. A choice of six atmospheric models (tropical, midlatitude 
summer and winter, subarctic summer and winter, and the U.S. 1962 
standard atmosphere) covering seasonal and latitudinal variations from 
sea level to 100 km are available, in addition to a capability of 
allowing the user to input atmospheric profile or meteorological para- 
meters of his own choosing. The program includes four aerosol models 
(average continental, urban, rural, and maritime) and permits either 
hazy (5-km visibility) or clear (23-km visibility) haze condition to be 
selected, in addition to the capability of selecting a particular 
visibility of the user's choosing. The model accounts for molecular 
absorption, molecular scattering, and aerosol extinction, plus atmos- 
pheric refraction and earth's curvature effects. 


LOWTRAN is presently the best available method for predicting 
atmospheric transmittance and, as such, should be incorporated in 
simulations and studies of electro-optic weapon/sensor systems per- 
formance. One advantage of LOWTRAN is that AFGL is continually improving 
and upgrading the computer code. For further documentation on LOWTRAN 
and its limitations, see Refs. 1-5, 11, and 14. 


» 


* BACKGROUND 


In examining AFGL's version of LOWTRAN (hereafter referred to 
AFGL/LOWTRAN), the author became aware that it would be possible to re- 
design the program without changing the basic methodology. This 
re-design achieves two gvals: 


a. a significant de’ ‘ease in program size and hence computer 
time; and 


b. a more structured program flow. 


The first goal ts critical for repeated uses of the program or if the 
program is to be incorporated into a weapon/sensor systems simulation. 
The second goal is tied to the first one in that a smooth, structured 
logic flow usually improves computer time. Moreover, a structured 


program is easier to understand, to use, and, when necessary, to modify 
and update. . 
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The author re-designed AFGL/LOWTRAN and produced the computer code 

‘ NWC/LOWTRAN. It should be emphasized that NWC/LOWTRAN does not attempt 
to modify the methodology or the algorithm used in AFGL/LOWTRAN in any 

way whatsoever. The only difference between the two versions is in the 

i computer coding of the algorithm. In fact, where possible, NWC/LOWTRAN 
uses the same variable names and statement labels as AFGL/LOWTRAN III 
(ref. 8). This was done to assist comparison between the two codes. 

| The author recently received the modifications to LOWTRAN III necessary 
to upgrade it to LOWTRAN IIIB status (Ref. 4). These changes have been 

implemented into NWC/LOWTRAN (Ref. 15). This report documents the 

| updated code, in addition to some corrections which have been brought to 

rt the author's attention (e.g., Refs. 9, 10, 14, and 16). 

{ 

} 


The author would particularly like to thank Mr. T. J. Smith of the 
Naval Weapons Center and Dr. A. Goroch of the Naval Environmental 


Prediction Research Facility for indicating several errors. In addition, 


the author would appreciate having any additional corrections to LOWTRAN 
brought to his attention at the following address: 


Commander (Code 3173) 
Attn: Dr. William M. Cornette 


i £ Naval Weapons Center 
| f: China Lake CA 93555 


Comm. Telephone (714) 939-3551 
Autovon 245-3551 
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, BASIC STRUCTURE 


The NWC/LOWTRAN computer code consists of four modules: 


(1) the main program; 


(2) the subroutine POINT, which computes the mean refractive i at. 
) |. index above and below a given altitude and interpolates ; 


exponentially to determine the bs oe ara absorber amounts 
at that altitude; 


(3) the subroutine ANGL, which calculates the initial zenith 
ey angle taking into account refraction effects; and 


eS NR AN ACER in BS HSA 


; . 
i | (4) a BLOCK DATA module, which inputs the spectral and 
transmittance data. 


, | 
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The number of executable statements contained in each module 


are: 
Main program 405 
. ; Subroutine POINT 26 
Subroutine ANGL 217 
BLOCK DATA _0 
TOTAL 648 


This is a reduction by over 200 lines of executable computer code 
from AFGL/LOWTRAN (Ref. 15). 


The actual run time for execution time of NWC/LOWTRAN is dependent 
upon the number of levels used in the geometry, whether the subroutine 
ANGL is called, and the width and resolution of the spectral band. 
However, the execution time for the sample output in Figure 1 (shown 
| later in this report) was approximately 470 milliseconds. 


INPUT /OUTPUT 


The basic inputs to NWC/LOWTRAN are only four cards: 


Card 1: MODEL, IAERO, IHAZE, JP, IM, Ml, M2, M3, ML, RO, VIS 
woen—~------ FORMAT (913, 3X,2F10. 3) 


MODEL selects one of the following model atmospheres: } 


. MODEL = O for a horizontal path when meteorological data is 
used. Instead of card 2, read altitude (km), pressure 
(mb), temperature (deg C), dew point temperature 
(deg C), relative humidity (%), water vapor density ! 


(gm. n°), ozone density (gn. an”), visibility (km), 
and range (km). P . 
ween en een-= FORMAT (3F10.3,2F5.2,2(1PE10. 3) ,2(OPF10.3)). 


MODEL = 1 specifies a tropical atmosphere. 

MODEL = 2 specifies a midlatitude summer atmosphere. 

MODEL = 3 specifies a midlatitude winter atmosphere. 

1 | MODEL = 4 specifies a sub-arctic summer atmosphere. 

|] MODEL = 5 specifies a sub-arctic winter atmosphere. . 
MODEL = 6 specifies the 1962 U.S. Standard atmosphere. - 

, | MODEL = 7 for a new model atmosphere (e.g. radiosonde data). 

a Read between cards 1 and 2, altitude (km), pressure (mb), 

temperature (deg C), dew point temperature (deg C), 


fbb ticietiot tcbenci t 2 


telative humidity (2%), water vapor density (gm. a), 
ozone density (gm. m~), and aerosol density (em™>), 


i TONE NRE HO 0g YET ee ae gir 
% rd BU Eg Ce 


me = . 
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ereewenennn= FORMAT (3F10.3,2FS.2,2(1PE10.31,2(OPE10.3)). 


Note that either dew point temperature, relative 
humidity, or water vapor density can be used. 


IAERO selects the type of aerosol attenuation: 


IAERO = 1 specifies an average continental aerosol model. 
IAERO = 2 specifies a rural aerosol model. 

IAERO = 3 specifies an urban aerosol model. 

IAERO = 4 specifies a maritime aerosol model. 


IHAZE selects the degree of aerosol attenuation: 


If IHAZE = 0 no aerosol scattering is computed. 

If IHAZE = 1 and VIS is non-zero, then aerosol attenuation 
for the visible range is used. 

If IHAZE = 1 or 2 and VIS is zero, then aerosol attenuation 
for 23 km and 5 km visible ranges, respectively, 
is used. 


JP selects the print option: 


JP = 0 for normal operation. 
JP = 1 to suppress the printing of the transmittance table. 


IM determines when radiosonde data is toa be read in: 
IM = 0 for normal operation or when subsequent calculations 
are to be run with MODEL = 7. 
IM = 1 for reading in radiosonde data initially for MODEL = 7. 


Ml, M2, and M3 are used to change temperature, H,0, and 0; altitude 


profiles, respectively, to another model value. 
ML is the number of levels of radiosonde data for MODEL = 7. 


RO is the input value for the radius of the earth. If RO = 0.0, 
then the program uses stored values. 


VIS is the visual range at sea level (km). 


Card 2: I = TYPE, LEN, H1, H2, ANGLE, RANGE, BETA---FORMAT (213,4X,5F10.3) 


ITYPE indicates the type of atmospheric path: 


ITYPE = 1, corresponds to a horizontal (constant pressure) 
path. Read Hl and RANGE. 
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ITYPE = 2, vertical or slant path between two altitudes. 
Read Hl and two other geometric parameters (e.g., H2 
and ANGLE). 

ITYPE = 3, vertical or slant path to space. 

* Read Hl and ANGLE. 


LEN selects the type of ray path to be used: 
LEN = 0 for normal operation of the program which selects the 
shorter path when applicable. 
LEN = 1 to select the longer path when applicable. 


observer altitude (km) 


H1 = 

H2 = source altitude (km) 

ANGLE = zenith angle at Hl (degrees) 

RANGE = path length (km) 

BETA = earth center angle (degrees) 

Card 3: V1, V2, DV--------------------- a ee ae em me FORMAT (3F10.3) 


Vl = initial frequency (wavenumber (cm-1)) value 

V2 = final frequency (wavenumber (cm-1)) value 

‘DV = frequency intervals at which transmittance is printed 
Note: Vl, V2, and DV must be integral multiples of 5 cm-l 


Card 4: IXY---------------------------=---------- FORMAT (13) 
IXY is the cycling indicator: 


° IXY = 0 to end data 
IXY = 1 for new card 3 only 
IXY = 2 to continue data 
IXY = 3 for new card 2 only 
IxY = 4 for new card 1 only 


A sample output is shown in Figure l. 
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PORTABILITY 


During the development of the NWC/LOWTRAN computer code, it was 
attempted to use a highly portable subset of ANSI FORTRAN (Ref. 12). To 
facilitate this, the PFORT verifier developed by Bell Telephone Labora- 


tories (Ref. 13) has been used to check the portability of the computer 
code. 


To the best of the author's knowledge, the only areas where diffi- 
culties may arise in using NWC/LOWTRAN with compilers and machines other 
than Univac compilers and computers, are as follows: 


(i) The program uses the Univac library routines ACOS (arccosine), 
ASIN (arcsine), and TAN (tangent); if not available on the desired 
machine (e.g., on Honeywell Series 32 computers) or if another name 
is used (e.g., ARCOS and ARSIN on the IBM 360 series computers), 

some changes will be necessary. 


(141) The values for the input/output units on the Univac 1110 
are 5 for the card reader and 6 for the line printer; these values 


are input to the program in the BLOCK DATA module for easy accessi- 
‘bility. 


(iii) Due to the large arrays used in NWC/LOWTRAN, much of the 
data in the BLOCK DATA module are assigned by DATA statements to an 
array name, or part of an array through an implicit DO; this 
feature may create certain difficulties on some compilers. 


If anyone experiences any problems using NWC/LOWTRAN due to com- 


piler or machine incompatibilities, the author would greatly appreciate 
being notified. 
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MAJOR SYMBOLS AND DEFINITIONS 


Absorption at frequency V; also average transmittance 

Aerosol number density 

Equivalent absorber amount per km at level J 

Angle of arrival at adjacent level 

Input zenith angle (degrees) 

Angle subtended at the earth's center as path traverses 
adjacent levels 

Total angle subtended by path at earth's center 

Equivalent absorber amount per km at level J + 1 

Conversion factor from degrees to radians 

Wavelength dependent coefficient used in refractive 
index expression 

Wavelength dependent coefficient used in refractive 
index expression 

Log absorption coefficient for water vapor 

Log absorption coefficient for uniformly mixed gases 

Log absorption coefficient for ozone 

Absorption coefficient for nitrogen (4m) 

Absorption coefficient for water vapor continuum 
(~4ym) 

Extinction coefficient for molecular: scattering 

Extinction coefficient for aerosol models 

Aerosol absorption coefficient 

Absorption coefficient for ozone (UV and visible regions) 

Water vapor amount (pr. cm/km) at level I 

Dew point temperature (°C) 

Path length from level I to Level I +1 

Wavenumber increment at which transmittance is calculated 

Height increment from level I to level I +1 : 

Equivalent absorber amounts per km at height Hl 

Equivalent absorber amount per km for water vapor at 
level Z(I) 

Equivalent absorber amount per km for carbon dioxide, 
etc. at level 2(I) 

Equivalent absorber amount per km for ozone at level Z(I) 

Equivalent absorber amount per km for nitrogen at level Z(I) 

Equivalent absorber amount per km for water vapor continuum 

? at level Z(1) 

Equivalent absorber amount per km for molecular scattering 
at level 2(1) 

Equivalent absorber amount per km for aerosol extinction 
at level Z(I) \ 

Equivalent absorber amount per km for ozone (UV and visible) © 
at level Z(I) 

Mean refractive index of layer above level Z(1) 
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Integrated absorber amount from level I to level I+1 

Factor for exponential and linear interpolation 

Transmission function logarithmic absorber amount scale 
for ozone 

Transmission function logarithmic absorber amount scale 
for water vapor and the uniformly mixed gases 

Altitude (km) 

Initial altitude (km) 


. Final altitude (km) 


Aerosol number density (no. cm_ 3) 
Estimated tangent height (km) 
Minimum altitude of path trajeetory (km) 


Aerosol number density (no. cn for 23 km visual range 


Aerosol number density (no. ca? ) for 5 km visual range 
Running integer used as altitude (level) indicator 
Indicator for type of aerosol model 

Number of levels in model atmosphere 


Frequency increment (em 4) 

Indicator for using subroutine ANGL 

Aerosol model indicator 

Parameter used when reading in a new atmospheric model 

Indicator for using subroutine POINT to calculate 
refractive index only (IP = 0) or equivalent ities 
amounts also (IP # 0). 

Card printer number 

Iteration counters 

Indicator for type of atmospheric path 

Frequency at which transmittance is calculated 

Starting frequency 

Last frequency 

Line printer number 

Parameter for terminating program and cycling indicator 

Running integer for altitude identification ~ 

Altitude indicator for minimum height of path 

Print option parameter 

Level indicator for ‘altitude Hl 


‘Level indicator for altitude H2 


Absorber indicator, K = 1, 2, 3, etc., corresponds to 
water vapor, uniformly mixed gases, ozone, etc., 
respectively 

Cycling parameter for downward looking paths 

Frequency indicator for ozone transmittance calculation 

Parameter used for defining longest of two paths 

Integer used to identify required model atmosphere 

Number of levels in radiosonde data input (MODEL - 7) 

Integer used to identify required model atmosphere 


— a 


P(M,I) 
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Integer for selecting temperature altitude profile for 
(M=M1) 

Integer for selecting water vapor altitude profile for 
(M=M2) 

Integer for selecting ozone altitude profile for (M=M3) 

Indicator for level below given input altitude used in 
POINT subroutine 

Number of levels in model atmosphere data 

Indicator for determining whether Hl or H2 coincide with 
leveis in the model atmosphere 

Value of NP for altitude Hl 

Value of NP for altitude H2 

Counters corresponding to WS1, WS2, WS3 

Pressure (mb) at level I for model atmosphere M 

Angle of arrival at H2 

3.141592654, that is 7 

Partial pressure of water vapor (in atmospheres) 

Total pressure in atmospheres 

Angular deviation of path from initial direction 

Product of total pressure (atm) and the square root of 
273/T(M,1) 

Path length (kn) 

Earth radius (kn) 

Refractive index of air at level I 

Relative humidity (%) 

Ratio of refractive indices of air above and below a 
given level 

Ratio of earth center distances between adjacent levels 

Earth radius (km) read in as input (=RE) 

The product of the sine of the initial zenith angle and 
the earth center distance to starting altitude 

Sine of angle of arrival at adjacent level 

Sine of the local zenith angle at a given level 

Slant range (km) 

Accumulated integrated absorption 

Temperature (K) for model atmosphere M at level I 

Zenith angle at a given level (in radians) 

Zenith angle at a given level (in degrees) 

Ambient temperature (°C) 

Transmittance scales for transmission functions 

Ratio of standard temperature (273.15K) to temperature 
at level I 2 

Ratio of 296.0K to temperature at — I % 

Ratio 273.15/(TMP + 273.15) 

Equivalent absorber amounts per km at a "given altitude 
obtained from POINT; also transmittance values at 
a given wavelength for each absorber type (K = 1,- 8) 

Total transmittance at frequency IV 

Absorption due to aerosol only at frequency IV 


V1 


v2 
W(K) 
WH(M, I) 


WL,WL1,WL2 
Wwo(M, 1) 


wsl 


WS2 


x 
XI 
XX 


x1 
X2 
YN 


YN1 
YN2 


YY 
Z(1),Z0(1) 
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Refractive index of layer above initial altitude Hl 

Refractive index of layer above final altitude H2 

Refractive index of layer above minimum altitude HMIN 

Integral of the equivalent absorber amounts from Hl to 
level I 

Visual range (km) at sea level 

Wavelength at which aerosol coefficients are read in (jim) 


Initial frequency for transmittance calculation, cat 


Final frequency for transmittance calculation, ca? 
Total equivalent absorber amount for entire path 
Water vapor density for atmospheric model M at level I 


-3 

(gm m ~) 
Wavelength in microns 
Ozone density for atmospheric model M at level I 
(gn n> 
Transmission function scaling factor for water vapor at 

given wavelength 
Transmission function scaling factor for carbon dioxide, 

etc., at given wavelength 
Transmission function scaling factor for ozone at given 

wavelength 
Input height to POINT subroutine 
Wavenumber interpolation parameter 
Wavenumber identification parameter for UV ozone trans- 

mittance calculation 
Earth center distance of level I 
Earth center distance of level I + l : 
Refractive index of layer below input height from POINT 

subroutine 
Refractive index of layer below initial altitude II1 
Refractive index of layer below final altitude H2 
Aerosol absorption coefficient at frequency V . ; 
Altitude at level I in km 
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CONCLUSIONS 


LOWTRAN is, in the author's opinion, the best available computer 
code for analyzing atmospheric transmittance effects on electro-optical 
systems. As such, it should be used in the evaluation of all electro- 
optic type sensors. With the release of LOWTRAN ‘IIIB, several previously 
valid criticisms have been accounted for; in addition, the U.S. Air 
Force Geophysics Laboratory continues to improve and upgrade the existing 
code. “4 


The author would appreciate being notified of any difficulties or 
problems that arise when using NWC/LOWTRAN. Also, any assistance that 
may be required, either in using the existing code or in modifying 
the code for use in a simulation, can be obtained by contacting the 
author. 
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